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The  optimum  design  of  thermoelectric  devices  (thermoelectric  generator  and  cooler)  in  connection  with 
heat  sinks  was  developed  using  dimensional  analysis.  New  dimensionless  groups  were  properly  defined 
to  represent  important  parameters  of  the  thermoelectric  devices.  Particularly,  use  of  the  convection  con¬ 
ductance  of  a  fluid  in  the  denominators  of  the  dimensionless  parameters  was  critically  important,  which 
leads  to  a  new  optimum  design.  This  allows  us  to  determine  either  the  optimal  number  of  thermocouples 
or  the  optimal  thermal  conductance  (the  geometric  ratio  of  footprint  of  leg  to  leg  length).  It  is  stated  from 
the  present  dimensional  analysis  that,  if  two  fluid  temperatures  on  the  heat  sinks  are  given,  an  optimum 
design  always  exists  and  can  be  found  with  the  feasible  mechanical  constraints.  The  optimum  design 
includes  the  optimum  parameters  such  as  efficiency,  power,  current,  geometry  or  number  of  thermocou¬ 
ples,  and  thermal  resistances  of  heat  sinks. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Thermoelectric  devices  (thermoelectric  generator  and  cooler) 
have  found  comprehensive  applications  in  solar  energy  conversion 
[1],  exhaust  energy  conversion  [2,3],  low  grade  waste  heat  recov¬ 
ery  [4-6],  power  plants  [7],  electronic  cooling  [8],  vehicle  air  con¬ 
ditioners,  and  refrigerators  [7],  The  most  common  refrigerant  used 
in  home  and  automobile  air  conditioners  is  R-134a,  which  does  not 
have  the  ozone-depleting  properties  of  Freon,  but  is  nevertheless  a 
terrible  greenhouse  gas  and  will  be  banned  in  the  near  future  [9], 
The  pertinent  candidate  for  the  replacement  would  be  thermoelec¬ 
tric  coolers.  Many  analyses,  optimizations,  even  manufacturers’ 
performance  curves  on  thermoelectric  devices  have  been  based 
on  the  constant  high  and  cold  junction  temperatures  of  the  devices. 
Practically,  the  thermoelectric  devices  must  work  with  heat  sinks 
(or  heat  exchangers).  It  is  then  very  difficult  to  have  the  constant 
junction  temperatures  unless  the  thermal  resistances  of  the  heat 
sinks  are  zero,  which  is,  of  course,  impossible. 

A  significant  amount  of  research  related  to  the  optimization  of 
thermoelectric  devices  in  conjunction  with  heat  sinks  has  been 
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conducted  as  found  in  the  literature  [10-30],  It  is  well  noted  from 
the  literature  that  there  is  the  existence  of  optimal  conditions  in 
power  output  or  efficiency  with  respect  to  the  external  load  resis¬ 
tance  for  a  thermoelectric  generator  (TEG)  or  the  electrical  current 
for  a  thermoelectric  cooler  (TEC).  Many  researchers  attempted  to 
combine  the  theoretical  thermoelectric  equations  and  the  heat  bal¬ 
ance  equations  of  heat  sinks,  and  then  to  optimize  design  parame¬ 
ters  such  as  the  geometry  of  heat  sinks  [10],  allocation  of  the  heat 
transfer  areas  of  heat  sinks  [12,13,18,19],  thermoelement  length 
[14],  the  number  of  thermocouples  [15],  the  geometric  ratio  of 
the  cross-sectional  area  of  thermoelement  to  the  length  [16],  and 
slenderness  ratio  (the  geometric  factor  ratio  of  n-type  to  that  of 
p-type  elements)  [17],  It  can  be  seen  from  the  above  literature  that 
the  geometric  optimization  of  thermoelectric  devices  is  important 
in  design  and  also  formidable  due  to  so  many  design  parameters. 
The  thermal  conductance  of  thermoelements  that  is  the  most 
important  geometric  parameter  has  been  often  addressed  in  anal¬ 
ysis,  which  is  the  product  of  three  parameters:  the  number  of  ther¬ 
mocouples,  the  geometric  ratio,  and  the  thermal  conductivity.  In 
order  to  reduce  the  optimum  design  parameters,  obviously  dimen¬ 
sionless  analyses  were  performed  in  the  literature  [21-26], 
Yamanashi  [21]  developed  optimum  design  introducing  dimen¬ 
sionless  parameters  for  a  thermoelectric  cooler  with  two  heat 
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Nomenclature 

A 

cross-sectional  area  of  thermoelement  (cm2) 

Tool 

temperature  of  fluid  2  (°C) 

A-, 

total  fin  surface  area  at  fluid  1  (cm2) 

fool.max 

maximum  temperature  of  fluid  1  (°C) 

a2 

total  fin  surface  area  at  fluid  2  (cm2) 

Tool  min 

minimum  temperature  of  fluid  1  (°C) 

Ab 

base  area  of  heat  sink  (cm2) 

Vn 

voltage  of  a  module  (V) 

COP 

the  coefficient  of  performance 

Wn 

power  output  (W) 

hr 

heat  transfer  coefficient  of  fluid  1  (W/m2  K) 

Wn 

power  input  (W) 

h2 

heat  transfer  coefficient  of  fluid  2  (W/m2  K) 

z 

the  figure  of  merit 

I 

electric  current  (A) 

L 

length  of  thermoelement  (mm) 

Creek  symbols 

k 

thermal  conductivity  (W/m  K),  k  =  kp  +  kn 

a 

Seebeck  coefficient  (V/K),  a  =  ap-ocn 

n 

the  number  of  thermocouples 

P 

electrical  resistivity  (£2  cm),  p  =  pp  +  pn 

Nk 

dimensionless  thermal  conductance  Nk  =  n(AklL)lt]2h2A2 

m 

fin  efficiency  of  heat  sink  1 

Nh 

dimensionless  convection,  Nb  =  rjthtA^ri2h2  A2 

f]2 

fin  efficiency  of  heat  sink  2 

N, 

dimensionless  current,  N/  =  a //(Ak/L) 

r]th 

thermal  efficiency  of  TEG 

Nv 

dimensionless  voltage,  Nv=  Vn/(naT^,2) 

Qi 

the  rate  of  heat  transfer  entering  into  TEG  (W) 

Subscripts 

02 

the  rate  of  heat  transfer  leaving  TEG  (W) 

P 

p-type  element 

Pa 

power  density  (W/cm2) 

n-type  element 

R 

electrical  resistance  of  a  thermocouple  (£2) 

opt 

optimal  quantity 

Rl 

load  resistance  of  a  thermocouple  (£2) 

1 /2opt 

half  optimal  quantity 

Rr 

dimensionless  resistance,  Rr  =  RJR 

h 

junction  temperature  at  fluid  1  (°C) 

Superscript 

t2 

junction  temperature  at  fluid  2  (°C) 

* 

dimensionless 

Tod 

temperature  of  fluid  1  (°C) 

sinks,  wherein  the  thermal  conductance  appears  twice  in  the  nom¬ 
inators  and  fourth  in  the  denominators  of  the  dimensionless 
parameters.  Although  his  work  led  to  a  new  approach  in  dimen¬ 
sionless  optimum  design,  the  analysis  encountered  difficulties  in 
optimizing  the  cooling  power  with  respect  to  the  thermal  conduc¬ 
tance  because  the  conductance  is  intricately  related  to  the  others. 
Later,  researchers  [22,24,25]  reported  optimum  design  using  the 
similar  dimensionless  parameters  used  by  Yamanashi  [21],  pre¬ 
senting  valuable  optimum  design  features  as  Xuan  [22]  optimized 
cooling  power  for  a  TEC  as  a  function  of  thermoelement  length,  Pan 
et  al.  [24]  showed  the  optimum  thermal  conductance  for  a  TEC 
with  a  given  cooling  power,  and  Casano  and  Piva  [25]  presented 
the  optimum  external  load  resistance  ratio  with  heat  sinks  for  a 
TEG  which  is  greater  than  unity.  There  are  also  some  experimental 
works  [27-30]  comparing  with  the  theoretical  thermoelectric 
equations.  Gou  et  al.  [27]  conducted  experiments  for  low-temper¬ 
ature  waste  heat  recovery  and  demonstrated  that  the  experimental 
results  were  in  fair  agreement  with  the  solution  formulas  origi¬ 
nally  derived  by  Chen  et  al.  [15]  from  the  general  theoretical 
thermoelectric  equations.  Chang  et  al.  [28]  and  Huang  et  al.  [29] 
conducted  experiments  for  a  TEC  from  a  heat  source  with 
air-cooling  and  water  cooling  heat  sinks,  respectively.  Casano  and 
Piva  [30]  reported  experimental  work  on  a  set  of  nine  thermoelec¬ 
tric  generator  modules  with  a  heat  source  on  one  side  and  a 
heat  sink  on  the  other  side.  After  deliberately  determined  the  heat 
leakage  which  turned  out  to  be  about  30%  of  the  supplied  heat 
source,  they  demonstrated  that  the  theoretical  performance  curves 
of  the  power  output  and  efficiency  as  a  function  of  the  external 
load  resistance  and  temperature  difference  were  in  good  agree¬ 
ment  with  the  measurements.  It  is  realized  from  the  above  exper¬ 
imental  works  that  the  theoretical  thermoelectric  equations  with 
the  heat  balance  equations  of  heat  sinks  can  reasonably  predict 
the  real  performance.  However,  proper  optimum  design  still 
remains  questionable. 

In  spite  of  many  efforts  for  optimum  design,  its  applications 
seem  greatly  challenging  to  system  designers  [1-3],  For  example, 
Hsu  et  al.  [3]  in  201 1  tested  an  exhaust  heat  recovery  system  both 


experimentally  using  an  automobile  and  mathematically  using 
computer  simulations.  They  found  a  reasonable  agreement  be¬ 
tween  the  measurements  and  the  simulations.  However,  they  ob¬ 
tained  the  power  output  of  12.41  W  over  24  thermoelectric 
generator  modules  with  the  exhaust  gas  temperature  of  573  K 
and  the  air  temperature  of  300  K.  When  the  power  output  was  di¬ 
vided  by  the  footprint  of  24  thermoelectric  generator  modules,  it 
gives  the  power  density  of  0.032  W/cm2,  which  seems  unusually 
small.  Karri  et  al.  [2]  in  2012  conducted  a  similar  experiment  with 
an  SUV  automobile.  This  time  they  designed  the  exhaust  heat 
recovery  system  with  an  optimum  coolant  flow  rate.  They  obtained 
the  power  output  of  550  W  over  16  thermoelectric  generator  mod¬ 
ules  with  the  exhaust  gas  temperature  of  686  K  and  the  coolant 
temperature  of  361  K,  which  provided  the  power  density  of 
0.61  W/cm2.  This  shows  a  significant  improvement,  indicating 
the  importance  of  optimum  design.  A  New  Energy  Development 
Organization  (NEDO)  Program  (Japan)  [7]  in  2003  also  reported  a 
similar  experiment  with  a  passenger  car,  obtaining  the  power  out¬ 
put  of  240  W  over  16  segmented-type  modules  with  the  exhaust 
gas  temperature  of  773  K  and  the  coolant  temperature  of  298  K, 
which  provided  the  power  density  of  ~1  W/cm2.  Notably,  the 
power  densities  obtained  are  no  way  to  evaluate  how  good  it  is  un¬ 
til  the  better  comes  because  proper  optimum  design  seems  not 
available. 

From  the  review  of  the  above  theoretical  and  experimental 
studies  including  optimum  design  in  the  literature,  it  is  summa¬ 
rized  that  the  proper  optimum  design  should  be  determined  basi¬ 
cally  not  only  by  the  power  output  for  TEG  (or  cooling  power  for 
TEC)  but  also  by  the  efficiency  (the  coefficient  of  performance) 
simultaneously  with  respect  to  both  the  external  load  resistance 
(or  the  electrical  current)  and  the  geometry  of  thermoelement 
which  refers  to  the  number  of  thermocouples  and  the  geometric 
ratio.  The  former  (external  load  resistance)  is  well  attained  in  the 
literature  but  the  latter  (geometry)  is  vague.  Therefore,  the  opti¬ 
mum  design  seems  incomplete.  This  is  the  rationale  why  the  pres¬ 
ent  paper  is  to  improve  the  optimum  design  introducing  new 
dimensionless  parameters. 
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Fluid  1 


1.1.  Thermoelectric  generator 

Let  us  consider  a  simplified  steady-state  heat  transfer  on  a  ther¬ 
moelectric  generator  module  (TEG)  with  two  heat  sinks  as  shown 
in  Fig.  la.  Each  heat  sink  faces  a  fluid  flow  at  temperature  Tm.  Sub¬ 
scripts  1  and  2  denote  hot  and  cold  quantities,  respectively.  We  as¬ 
sume  that  the  electrical  and  thermal  contact  resistances  in  the  TEG 
are  negligible,  the  material  properties  are  independent  of  temper¬ 
ature,  and  also  the  TEG  is  perfectly  insulated.  The  TEG  has  a  num¬ 
ber  of  thermocouples,  of  which  each  thermocouple  consists  of  p- 
type  and  n-type  thermoelements  with  the  same  dimensions  as 
shown  in  Fig.  lb.  It  is  noted  that  the  thermal  resistance  of  heat  sink 
1  can  be  expressed  by  the  reciprocal  of  the  convection  conductance 
»7ihiAi,  where  r\ i  is  the  fin  efficiency,  h,  is  the  convection  coeffi¬ 
cient,  and  is  the  total  surface  area  in  the  heat  sink  1.  We  hereaf¬ 
ter  use  the  convection  conductance  rather  than  the  thermal 
resistance.  The  basic  equations  for  the  TEG  with  two  heat  sinks 
are  given  by 


_  n(Ak/L ) 


(6) 


The  dimensionless  convection  is  an  important  geometry  of  heat 
sinks.  Since  the  convection  conductance  at  fluid  2  is  given  as  men¬ 
tioned  before,  the  dimensionless  convection,  the  ratio  of  convec¬ 
tion  conductance  in  fluid  1  to  fluid  2,  is  defined  by 


(7) 


The  dimensionless  electrical  resistance,  the  ratio  of  the  load 
resistance  to  the  electrical  resistance  of  thermocouple,  is  given  by 


(8) 


Since  the  fluid  temperature  T„,2  at  fluid  2  is  given  as  mentioned 
before,  the  dimensionless  temperatures  are  defined  by 


Qi 

^MiCU-T,) 

(1) 

Qi 

=  n(u '.IT,-ll2R  +  ~(h 

A 

(2) 

q2 

=  n(aIT2  +  ^I2R  +  ^(h 

-">) 

(3) 

0.2 

=  ri2h2A2(T2  -  T„ 2) 

(4) 

1  = 

a  (7,  -  T2) 

Rl  +  R 

(5) 

where  a  =  <xp  -  a„,  k  =  kp  +  kn,  and  p  =  pp  +  pn.  Eqs.  (1 )— (5)  can  be 
solved  for  Ti  and  T2,  providing  the  power  output.  However,  in  order 
to  study  the  optimization  of  the  TEG,  several  dimensionless  param¬ 
eters  are  introduced.  As  mentioned  in  Section  of  Introduction,  it  is 
reminded  that  optimum  design  should  consider  not  only  the  power 
output  but  also  the  efficiency  simultaneously  with  respect  to  both 
the  external  load  resistance  and  the  geometry  of  thermoelement 
which  refers  to  the  number  of  thermocouples  and  the  geometric  ra¬ 
tio.  In  order  to  reveal  the  effect  of  thermal  conductance  n(AklL),  the 
thermal  conductance  is  placed  in  the  nominator  while  the  convec¬ 
tion  conductance  rj2h2A2  in  fluid  2  is  placed  in  the  denominator  of 
the  parameter.  The  convection  conductance  r)2h2A2  and  tempera¬ 
ture  7^2  at  fluid  2  are  assumed  to  be  given.  The  dimensionless  ther¬ 
mal  conductance,  the  ratio  of  thermal  conductance  to  the 
convection  conductance  in  fluid  2,  is  defined  by 


7, 

7^ 

(9) 

A 

(10) 

7ool 

=  7^2 

(11) 

The  dimensionless  power  and  heat  transfer  are  defined  by 
dividing  by  the  product  of  the  convection  conductance  and  the 
temperature  of  fluid  2  so  that  the  quantities  depend  only  on  the 
nominators  not  on  the  denominators  since  the  denominator  is  as¬ 
sumed  to  be  constant  or  given.  The  two  dimensionless  rates  of  heat 
transfer  and  the  dimensionless  power  output  are  defined  by 


q;  m 

w:  = 


Q 1 

Q2 

Wn 


(12) 

(13) 

(14) 


It  is  noted  that  the  above  dimensionless  parameters  are 
based  on  the  convection  conductance  in  fluid  2,  which  means  that 
rj2h2A2  Too2  should  be  initially  provided.  Also  note  that  the 
thermal  conductance  n(Ak/L)  appears  only  in  Eq.  (6)  among  other 
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parameters  so  that  the  thermal  conductance  can  be  examined  for 
optimization.  Using  the  dimensionless  parameters  defined  in 
Eqs.  (6)-(ll),  Eqs.  (l)-(5)  reduce  to  two  formulas  as: 


*(r-r,)_zr^(r,-r1)r,_zr^(rt-n) 

Nk  Rr  +  1  2(Rr  + 1)2  V1  2)1 

t*  - 1  zt,x;2 (t;  -  r2)T2  zt,x2 (t;  -  r*): 

Nk  '  2(Rr  +  l)2 


where  Z  is  called  the  figure  of  merit  (Z  =  a2/pk).  Eqs.  (15)  and  (16) 
can  be  solved  for  T\  and  V2.  The  dimensionless  temperatures  are 
then  a  function  of  five  independent  dimensionless  parameters  as 


TJ  =f(Nk,Nh,Rr,roa,Zroo2)  (17) 

r2  =f(Nk,Nh,RnrDO,zroo2)  (is) 

Tx,  is  the  input  and  ZT„,2  is  the  material  property  with  the  input, 
and  both  are  initially  provided.  Therefore,  the  optimization  can  be 
performed  only  with  the  first  three  parameters  ( Nk ,  Nh}  and  Rr). 
Once  the  two  dimensionless  temperatures  (T)  and  V2)  are  solved 
for,  the  dimensionless  rates  of  heat  transfer  at  both  hot  and  cold 
junctions  of  the  TEG  can  be  obtained  as: 


Q.\=Nh{roo-r, )  (19) 

0.2  =  T2  -  1  (20) 

Then,  we  have  the  dimensionless  power  output  as 

w;  =  q;-02  (21) 

Accordingly,  the  thermal  efficiency  is  obtained  by 


Defining  Nj  =  txIL/Ak,  the  dimensionless  current  is  obtained  by 


Also,  defining  Nv=Vlnu,Too2 
tained  by 


Nv  = 


N,Nk 


the  dimensionless  voltage  is  ob- 


(24) 


With  the  inputs  (T^  and  ZT^2),  we  begin  developing  the  optimi¬ 
zation  with  the  dimensionless  parameters  (Nfc,  Nh,  and  Rr)  itera¬ 
tively  until  they  converge.  It  is  found  that  both  Nk  and  Rr  show 
their  optimal  values  for  the  dimensionless  power  output,  while 
Nh  does  not  show  the  optimal  value  showing  that  the  dimension¬ 
less  power  output  monotonically  increases  with  increasing  Nh.  This 
implies  that,  if  Nh  is  given,  the  optimal  combination  of  Nk  and  Rr 
can  be  obtained.  However,  the  dimensionless  convection  Nh  actu¬ 
ally  presents  the  feasible  mechanical  constraints.  Thus,  we  first 
proceed  with  a  typical  value  of  Nh  =  1  for  illustration  and  later 
examine  the  variety  of  Nh  with  some  practical  design  examples. 

Suppose  that  we  have  two  initial  inputs  of  T^  =  2.6  (two  fluid 
temperatures)  and  ZT^2  =  1-0  (materials)  along  with  Nh=  1.  We 
then  determine  the  optimal  combination  for  Nk  and  Rr,  which 
may  be  obtained  either  graphically  or  using  a  computer  program. 
We  first  use  the  graphical  method  at  this  moment  and  later  the 
program  for  multiple  computations.  The  dimensionless  power  out¬ 
put  W*  and  thermal  efficiency  rjth  are  together  plotted  as  a  function 
of  Rr,  which  are  presented  in  Fig.  2a.  Both  W*  and  rjth  with  respect 
to  Rr  indeed  show  their  optimal  values  that  appear  close.  We  are 
interested  primarily  in  the  power  output  and  secondly  in  the  effi¬ 
ciency.  However,  since  they  are  close  each  other,  we  herein  use  the 
power  output  for  the  optimization.  It  should  be  noted  that  the 
dimensionless  maximum  power  output  does  not  occur  at  Rr  =  1 


from  Fig.  2a  as  usually  assumed  for  a  TEG  without  heat  sinks,  but 
approximately  at  Rr  =  1.7,  because  the  dimensionless  temperatures 
T‘t  and  V2  in  Fig.  2b  are  no  longer  constant.  This  is  often  a  confusing 
factor  in  optimum  design  with  a  TEG  with  two  heat  sinks.  We 
should  not  assume  that  Rr  is  equal  to  unity  for  a  TEG  with  heat 
sinks. 

With  the  dimensionless  parameters  obtained  (Nh  =  l,  Rr=  1.7, 
T^  =  2.6,  and  ZT^  =  1.0),  we  now  plot  the  dimensionless  power 
output  W*  as  a  function  of  the  dimensionless  thermal  conductance 
Nk  defined  in  Eq.  (6)  along  with  the  thermal  efficiency  rjth ,  which  is 
shown  in  Fig.  3a.  We  find  an  optimum  W*  approximately  at 
Nk  =  0.3.  Actually,  the  optimal  values  of  Nk  and  Rr  should  be  iterated 
until  the  two  simultaneously  converge.  From  Nk  =  n{AklL)lr\2h2A2  as 
shown  in  Eq.  (6),  the  Nk  actually  determines  the  number  of  thermo¬ 
couples  n  if  the  geometric  ratio  A/L  and  r]2h2A2  are  given  or  vice 
versa.  The  dimensionless  power  output  W*  first  increases  and  later 
decreases  with  increasing  Nk.  It  is  important  to  realize  that,  if  rj2h2 
A2  is  given,  there  is  an  optimal  number  n  of  thermocouples  (or 
optimal  thermal  conductance  AkjL)  in  the  thermoelectric  module, 
which  is  usually  unknown.  Physically,  the  surplus  number  of  ther¬ 
mocouples  virtually  increases  the  thermal  conduction  more  than 
the  production  of  power,  which  causes  the  net  power  output  to  de¬ 
cline.  There  is  another  important  aspect  of  the  optimal  dimension¬ 
less  thermal  conductance  of  Nk  =  0.3,  which  is  that  the  module 
thermal  conductance  riAkIL  directly  depends  on  the  r\2h2A2.  In 
other  words,  the  module  thermal  conductance  nAk/L  must  be 
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Fig.  3.  (a)  Dimensionless  power  output  W’n  and  thermal  efficiency  ijtfl  versus 
dimensionless  thermal  conductance  Nk  and  (b)  high  and  low  junction  temperatures 
(F  and  T'? ;  versus  dimensionless  thermal  conductance  Nk.  These  plots  were 
generated  using  Nh  =  1,  Rr=  1.7,  F  =  2.6  and  ZT,2  =  1.0. 


Nh 


Fig.  4.  Optimal  dimensionless  power  output  W’opt  and  efficiency  popt  versus 
dimensionless  convection  Nh.  This  plot  was  generated  using  F  =  2.6  and 

2Tqo2  =  1 .0. 


1l2h2A2  (W/K) 


Fig.  5.  Optimal  power  output  Wopt  versus  convection  conductance  ;/ 2/12A2  in  fluid  2 
as  a  function  of  dimensionless  convection  Nh.  This  plot  was  generated  using 
Tx2  =  25  °C,  F  =  2.6  and  ZJ„_2  =  1.0. 


redesigned  on  the  basis  of  the  r]2h2A2  to  meet  Nk  =  0.3.  The  infor¬ 
mation  of  the  optimal  thermal  conductance  [(n)(A/I)(k)]  is  particu¬ 
larly  important  in  design  of  microstructured  or  thin-film 
thermoelectric  devices.  Furthermore,  there  is  a  potential  to  im¬ 
prove  the  performance  or  to  provide  the  variety  of  the  geometry 
by  reducing  the  thermal  conductivity  fc. 

The  dimensionless  high  and  low  junction  temperatures  are  pre¬ 
sented  in  Fig.  3b.  As  Nk  decreases  towards  zero,  T]  and  V2  approach 
Tx  and  1,  respectively.  This  indicates  that  the  thermal  resistances 
of  two  heat  sinks  approaches  zero,  which  never  happens.  It  is  noted 
that  the  thermal  efficiency  approaches  the  theoretical  maximum 
efficiency  of  0.2  for  the  given  fluid  temperatures  as  Nk  approaches 
zero. 

Since  there  is  an  optimal  combination  of  Nk  and  Rr  for  a  given 
Nh,  we  can  plot  the  optimal  dimensionless  power  output  W*pt 
and  optimal  thermal  efficiency  t]opt  as  a  function  of  dimensionless 
convection  Nk,  which  is  shown  in  Fig.  4.  It  is  very  interesting  to 
note  in  Fig.  4  that,  with  increasing  Nh,  r\opt  barely  changes,  while 
W*opt  monotonically  increases.  According  to  Eq.  (14),  the  actual 
optimal  power  output  Wopt  is  the  product  of  W*opl  and  ij2h2A2, 
seemingly  increasing  linearly  with  i)Jd2A2-  In  practice,  there  is  a 
controversial  tendency  that  Nh  may  decrease  systematically  with 
increasing  ri2h2A2  if  h,A,  is  limited.  As  a  result  of  this,  it  is 
needed  to  examine  the  variety  of  the  ij2h2A2  as  a  function  of  Nh 
for  the  optimal  power  output.  Fig.  5  reveals  the  intricate  relation¬ 
ship  between  rj2h2A2  and  Nh  (or  t/1hiA1)  along  with  the  optimum 


actual  power  output  (not  dimensionless),  which  would  lead  system 
designers  to  a  variety  of  possible  allocations  (772/12  A2  and  Nh)  for 
their  optimal  design.  Now  we  look  into  the  actual  optimal  design 
with  the  actual  values. 

For  example,  using  Figs.  4  and  5,  we  develop  an  optimal  design 
for  automobile  exhaust  gas  waste  heat  recovery.  A  thermoelectric 
generator  module  with  a  5-cm  x  5-cm  base  area  is  subject  to  ex¬ 
haust  gases  at  500  °C  in  fluid  1  and  air  at  25  °C  in  fluid  2.  We  esti¬ 
mate  an  available  maximum  convection  conductance  in  fluid  1 
(exhaust  gas)  with  771=0.8,  h,  =  60  W/m2  K,  and  A,  =  1000  cm2 
and  also  an  available  maximum  convection  conductance  in  fluid 
2  (air)  with  rj2  =  0.8,  h2  =  60  W/m2  K,  and  A2  =  1000  cm2,  which 
gives  77ihiAi  =  772/I2A2  =  4.8  W/K  and  Nh=  1.  Note  that  >71  and  r\2 
are  typical  fin  efficiencies,  and  and  h2  are  the  reasonable  convec¬ 
tion  coefficients  for  exhaust  gas  heating  and  air  cooling,  of  which 
the  convection  coefficients  with  exhaust  gas  or  air  typically  have 
values  ranging  between  20  and  100  W/m2  K  depending  on  the  flow 
rate  and  the  type  of  fin.  The  middle  value  of  60  W/m2  K  was  used  in 
the  present  work.  Each  area  of  A,  and  A2  is  based  on  20  fins  (two 
sides)  with  a  fin  height  of  5  cm  on  a  5-cm  x  5-cm  base  area  of 
the  module,  which  gives  an  total  fin  area  (5  cm  x  5  cm  x  2 
sides  x  20  fins  =  1000  cm2).  The  typical  thermoelectric  material 
properties  are  assumed  to  be  ap  =  -an  =  220  pV/K,  pp  =  p„  =  1.0 
x  10  3  £1  cm,  and  kp  =  k„  =  1.4  x  10  2  W/cm  K.  The  above  data 
approximately  determines  three  dimensionless  parameters  as 
Nh  =  1,7^  =  2.6  and  ZT ^  =  1.0. 
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Inputs  and  results  from  the  dimensional  analysis  for  a  TEG. 


Inputs 

Dimensionless  (w*  opt) 

Actual  (WnJ,pt) 

Too,  =  500  '  C,  T„i2  =  25  °C,  AT„  =  475  °C 

IV*  =  0.3 

n  =  254 

A  =  2  mm2, 1  =  1  mm 

N„=  1 

f/,1 h  A,  =  4.8  W/K 

J/2  =  0.8,  h2  =  60  W/m2  K,  A2  =  1000  cm2 

Rr  =  1.7 

RL  =  1.7  x  n  x  R  =  4.32  Q 

rt2h2A2  =  4.8  W/K 

C  =  2.6 

Toot  =  500  °C 

Base  area  Ab  of  module  =  5  cm  x  5  cm 

ZF,;2  =  1.0 

ZT002  =  1.0 

ap  =  — a„  =  220  pV/K 

r;  =  2.172 

T,  =  374  °C 

Pp  =  Pn  =  1-0  x  10-3ncm 

T‘2  =  1.367 

T2  =  137  °C 

kp  =  k„  =  1.4  x  10-2  W/cm  K 

WJ  =  0.045 

W„  =  65.0  W 

(Z  =  3.457  x  fijp’lC1) 

rfrh  =  0.108 

nth  =  0.108 

(J?  =  0.01  Cl  per  thermocouple) 

N,  =  0.306 

1  =  3.9  A 

(»/ 1  =  0.8,  h,  =  60  W/m2  K,  A,  =  1000  cm2) 

Nv  =  0.50 

V=  16.7  V 

(Power  density  P„  =  W„IAb) 

- 

Pd  =  2.6  W/cm2 

As  mentioned  before,  the  present  dimensional  analysis  enables 
the  three  dimensionless  parameters  (Nh=  1,  —  2.6  and 

ZToo 2  =  1.0)  to  determine  the  rest  two  optimal  parameters,  which 
are  found  to  be  Nk  =  0.3  and  Rr=  1.7  as  shown  before.  This  leads 
to  a  statement  that,  if  two  individual  fluid  temperatures  on  heat 
sinks  connected  to  a  thermoelectric  generator  module  are  given, 
an  optimum  design  always  exists  with  the  feasible  mechanical 
constraints  that  present  N/,.  This  optimal  design  is  indicated 
approximately  at  Point  1  in  Fig.  5.  Note  that  there  are  ways  to  im¬ 
prove  the  optimal  power  output,  increasing  either  r]2h2A2  or  Nh  or 
both,  which  obviously  depends  on  the  feasible  mechanical  con¬ 
straints,  whichever  is  available. 

The  inputs  and  optimum  results  at  Point  1  in  Fig.  5  are  summa¬ 
rized  in  Table  1.  The  inputs  are  the  geometry  of  thermocouple,  the 
material  properties,  two  fluid  temperatures,  and  the  available  con¬ 
vection  conductance  in  fluid  2.  The  dimensionless  results  are  con¬ 
verted  to  the  actual  quantities  as  shown.  The  maximum  power 
output  is  found  to  be  65.0  W  for  the  5  cm  x  5  cm  base  area  of  the 
module.  The  power  density  is  calculated  to  be  2.6  W/cm2,  which 
appears  significantly  high  compared  to  an  available  power  density 
of  ~  1  W/cm2  with  the  similar  operating  conditions  by  NEDO  pro¬ 
gram  (Japan)  [7], 

Air  was  used  in  fluid  2  so  far.  However,  we  want  to  see  the  effect 
of  r]2h2A2  or  Nh  by  changing  fluid  2  from  air  to  liquid  coolant. 
Otherwise  the  same  conditions  were  applied  to  as  the  previous 
example.  We  then  estimate  an  available  convection  conductance 
in  fluid  1  (exhaust  gas)  with  the  same  one  of  rj i  =  0.8,  h i  =  60 
W/m2  K,  and  At  =  1000  cm2,  but  an  available  convection  conductance 
in  fluid  2  (liquid  coolant)  with  172  =  0.8,  h2  =  3000  W/m2  K,  and 
A2  =  100  cm2,  which  gives  i/jhjAt  =  4.8  W/K  and  r\2h2A2  =  24  W/K, 
respectively,  which  yields  Nh  =  0.1.  The  area  of  A,  is  based  on  20 
fins  (two  sides)  with  a  fin  height  of  5  cm  for  the  5-cm  x  5-cm  base 
area  of  the  module  and  A2  is  estimated  to  be  one  tenth  of  A,  (liquid 
coolant  does  not  require  a  large  heat  transfer  area).  These  inputs 
and  optimum  results  give  all  the  five  dimensionless  parameters 
as  Nk  =  0.07,  Nh  =  0.1,  Rr  =  1.5,  T^  =  2.6  and  ZT^  1.0,  for 
which  the  optimum  at  r)2h2A2  =  24  W/K  is  indicated  at  Point  2  in 
Fig.  5.  The  effect  of  Nh  on  the  high  and  cold  junction  temperatures 
was  also  presented  in  Fig.  6.  It  is  interesting  to  see  that,  although  a 
small  variation  in  the  optimal  power  outputs  between  Point  1 
(N/,  =  1)  and  Point  2  (Nh  =  0.1)  appears  in  Fig.  5,  a  significant  tem¬ 
perature  variation  between  Nh  =  1  and  Nh  =  0.1  appears  in  Fig.  6. 
This  may  be  an  important  factor  particularly  when  thermoelectric 
materials  are  considered  in  the  optimal  design.  The  proximity  of 
the  power  outputs  between  Points  1  and  2  is  an  example  showing 
the  variety  of  the  mechanical  constraints  (i/ihiAi  and  r]2h2A2)  even 
with  the  same  power  outputs.  It  is  important  to  realize  that,  when 
»71h1A1  is  limited,  simply  increasing  i]2h2A2  invokes  decreasing  Nh, 
which  results  in  decreasing  not  only  the  high  and  cold  junction 
temperatures  but  also  slightly  the  temperature  difference  as 


Nh 


Fig.  6.  Hot  and  cold  junction  temperatures  and  optimal  efficiency  versus  dimen¬ 
sionless  convection.  This  plot  was  generated  with  T„2  =  25  °C,  F  =  2.6  and 

zr„2  =  1.0. 

shown  in  Fig.  6.  The  coexistence  that  the  Seebeck  coefficient  de¬ 
creases  with  decreasing  the  temperature  and  reducing  the  temper¬ 
ature  difference  diminishes  the  performance  will  cause  the  power 
output  to  decline.  However,  increasing  t]2h2A2  will  directly  increase 
the  power  output  as  mentioned  earlier.  The  net  power  output  of 
loss  and  gain  by  increasing  ri2h2A2  may  be  a  role  of  system  de¬ 
signer.  Anyhow  there  will  be  a  small  change  in  the  efficiency. 

2.  Thermoelectric  cooler 

Let  us  consider  a  simplified  steady-state  heat  transfer  on  a  ther¬ 
moelectric  cooler  module  (TEC)  with  two  heat  sinks  as  shown  in 
Fig.  7.  Each  heat  sink  faces  a  fluid  flow  at  temperature  T Sub¬ 
scripts  1  and  2  denote  the  entities  of  fluids  1  and  2,  respectively. 
Consider  that  an  electric  current  is  directed  in  a  way  that  the  cool¬ 
ing  power  Ch  enters  heat  sink  1.  We  assume  that  the  electrical  and 
thermal  contact  resistances  in  the  TEC  are  negligible,  the  material 
properties  are  independent  of  temperature,  and  also  the  TEC  is  per¬ 
fectly  insulated.  The  TEC  has  a  number  of  thermocouples,  of  which 
each  thermocouple  consists  of  p-type  and  n-type  thermoelements 
with  the  same  dimensions. 

The  basic  equations  for  the  TEC  with  two  heat  sinks  are  given  by 


a,  =J7,Mi(r0ol-r1) 

(25) 

a1=n^a/T1-jf2R  +  ^(T,-r2)) 

(26) 

Q2  =  n(a/T2  +  l?R+%  -  r2)) 

(27) 

Q2  =  r)2h2A2(T2  -  Tool) 

(28) 
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Fluid  1 


Fig.  7.  Thermoelectric  cooler  module  (TEC). 


where  a  =  ap  -  <x„,  k  =  kp  +  kn,  and  p  =  pp  +  pn.  In  order  to  study  the 
optimization  of  the  TEC,  several  dimensionless  parameters  are 
introduced.  The  dimensionless  thermal  conductance,  which  is  the 
ratio  of  thermal  conductance  to  the  convection  conductance  in  fluid 
2,  is 


mm 

N  ~  JfeM2 

(29) 

The  dimensionless  convection,  which  is  the  ratio  of  convection  con- 

ductance  in  fluid  1  to  fluid  2,  is 

=  |M, 
h  r,2h2A2 

(30) 

The  dimensionless  current  is  given  by 

(31) 

The  dimensionless  temperatures  are  defined  by 

T'-b 

(32) 

(33) 

r"  =  fe 

(34) 

The  dimensionless  cooling  power,  rate  of  heat 
electrical  power  input  are  defined  by 

liberated  and 

tv-  Qi 

1  n2h2A2Tx2 

(35) 

<r- 

i\2h2A2T.y2 

(36) 

W„ 

w*  -  n 

"  tl2h2A2T,x2 

(37) 

It  is  noted  that  the  above  dimensionless  parameters  are  based 
on  the  convection  conductance  in  fluid  2,  which  means  that  rj2h2. 
A2Too2  should  be  initially  provided.  Using  the  dimensionless  param¬ 
eters  defined  in  Eqs.  (29)-(34),  Eqs.  (25)-(28)  reduce  to  two 
formulas  as: 


NhK-T\) 

Nk 


=  n,t ;  - 


N f 

2ZJx2 


+  C n-r2) 


TJ  -  •  1 

% 


=  N,r2  + 


N jf 

2 ZJX 


+  iri  ~  t2) 


(38) 

(39) 


Eqs.  (38)  and  (39)  can  be  solved  for  T]  and  V2.  The  dimensionless 
temperatures  are  then  a  function  of  five  independent  dimensionless 
parameters  as 

T*  =f(Nk.Nh,Nhrx,ZJoo2)  (40) 

r2=f(Nk,Nh,NIJl,ZJoo2)  (41) 

T^  is  the  input  and  ZTx2  is  the  material  property  with  the  input, 
and  both  are  initially  provided.  Therefore,  the  optimization  can  be 
performed  only  with  the  first  three  parameters  (Nk,  Nj,,  and  N;). 
Once  the  two  dimensionless  temperatures  (7j  and  T‘)  are  solved 
for,  the  dimensionless  rates  of  heat  transfer  at  both  junctions  of 
the  TEC  can  be  obtained  as: 


Q-]  =  Nh{Tmx-ri)  (42) 

Q*  =  T*  -  1  (43) 

Qj  is  called  the  dimensionless  cooling  power.  Then,  we  have  the 
dimensionless  power  input  as 

W*„  =  Q*  -  Q;  (44) 

Accordingly,  the  coefficient  of  performance  is  obtained  by 


(45) 


Defining  Nv  =  the  dimensionless  voltage  is  obtained  by 


Nv 


W*n 

N,Nk 


(46) 


With  the  inputs  (T^  and  ZJcxi2),  we  try  to  find  the  optimal  com¬ 
bination  for  the  dimensionless  parameters  ( Nk,Nh ,  and  /V,)  itera¬ 
tively  until  they  converge.  It  is  found  that  both  Nk  and  Nj  show 
the  optimal  values  for  the  dimensionless  cooling  power  Q) ,  while 
Nh  does  not  show  the  optimal  value  showing  that  the  dimension¬ 
less  cooling  power  Qj  monotonically  increases  with  increasing 
Nh.  This  implies  that,  if  any  Nh  is  given,  the  optimal  combination 
of  Nk  and  JV;  can  be  obtained.  However,  the  dimensionless  convec¬ 
tion  Nh  actually  presents  the  feasible  mechanical  constraints.  Thus, 
we  proceed  with  a  typical  value  of  Nh  =  1  for  illustration  and  later 
examine  the  variety  of  Nh  with  a  practical  design  example. 

Suppose  that  we  have  T,%  =  0.967  (two  arbitrary  fluid  tempera¬ 
tures)  and  ZTa o2  =  1.0  (materials)  along  with  Nh  =  1  as  inputs.  Then, 
we  can  determine  the  optimal  combination  for  N,  and  Nk,  which 
may  be  obtained  either  graphically  or  using  a  computer  program. 
We  first  use  the  graphical  method  at  this  moment  and  later  the 
program  for  multiple  computations  (a  Mathematical  software 
Mathcad  was  used).  The  optimal  combination  of  N(  and  Nk  for  each 


Fig.  8.  Dimensionless  cooling  power  Qj,  power  input  W'n  and  COP  versus 
dimensionless  current  N,.  This  plot  was  generated  with  Nt  =  0.3,  Nh  =  1, 
=  o.967  and  zr«a= 1.0. 
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Fig.  9.  Dimensionless  cooling  power  Qj  and  COP  versus  dimensionless  thermal 
conductance  Nk.  This  plot  was  generated  with  Nk  =  1,  N,  =  0.5,  T^  =  0.967  and 


maximum  cooling  power  are  found  to  be  N,  =  0.5  and  Nk=  0.3, 
respectively,  which  are  shown  in  Figs.  8  and  9.  The  maximum  cool¬ 
ing  power  of  Qi  =  0.037  in  both  figures  is  actually  the  optimal 
dimensionless  cooling  power  Qj  opt.  However,  the  COP  also  shows 
an  optimal  value  at  N(  =  0.074,  which  gives  Qj  =  0.006.  The  optimal 
COP  usually  gives  a  very  small  cooling  power  or  sometimes  even  no 
exists,  which  seems  impractical,  albeit  the  high  COP.  Therefore,  it  is 
needed  to  have  a  practical  point  for  the  optimal  COP,  which  is 
determined  in  the  present  work  to  be  the  midpoints  of  the  opti- 
mumN;  and  Nk.  For  example,  the  practical  optimal  COP  in  this  case 
occurs  simultaneously  at  N,=  0.25  and  Nk  =  0.15,  which  leads  to 
Qj  =  0.019  that  may  be  seen  after  re-plotting  with  the  two  values. 

The  existence  of  an  optimum  cooling  power  as  a  function  of  cur¬ 
rent  is  a  well  known  characteristic  of  TECs.  However,  the  existence 
of  the  optimum  Nk  in  TECs  has  not  been  found  in  the  literature  to 
the  author’s  knowledge.  With  Eq.  (29)  that  is  Nk  =  n(Ak/L)lri2h2A2, 
the  optimum  of  Nk  =  0.3  implies  that  the  module  thermal  conduc¬ 
tance  n(AklL)  is  at  optimum  since  the  t]2h2A2  is  given,  which  leads 
to  the  optimum  n  (the  number  of  thermocouples)  if  Ak/L  is  given  or 
vice  versa.  This  is  one  of  the  most  important  optimum  processes  in  de¬ 
sign  of  a  thermoelectric  cooler  module.  It  is  good  to  know  in  Fig.  10 
that  the  dimensionless  temperature  T\  becomes  lowest  at  the  optimal 
dimensionless  cooling  power  Qj ,  not  at  the  optimum  COP. 

We  also  consider  two  heat  sinks  as  a  unit  without  a  TEC  to 
examine  the  limitation  of  use  of  the  TEC,  which  is  shown  in 
Fig.  11.  The  geometry  of  the  unit  is  the  same  as  the  one  shown  in 
Fig.  7  except  that  there  is  no  TEC  between  the  heat  sinks.  There 
should  be  a  cooling  rate  with  given  fluid  temperatures,  which  is 


Nk 


Fig.  10.  Dimensionless  temperatures  versus  dimensionless  thermal  conductance. 
This  plot  was  generated  with  Nk  =  1,  N,  =  0.5,  T^  =  0.967  and  =  1.0. 


Qo.  We  want  to  compare  the  cooling  power  Q,  with  this  cooling 
rate  Qo  to  determine  the  limit  of  use  of  the  TEC. 

The  basic  equations  for  the  unit  can  be  expressed  as 


Qjo  =  hi  Ax  (Txr  -  T0)  (47) 

Q.0  =  ihhiAjiJa  —  Tata)  (48) 

The  dimensionless  groups  for  the  unit  are 


T°  =  ft2  (49) 

<50> 

Using  Eqs.  (30),  (34),  (49)  and  (50),  Eqs.  (47)  and  (48)  reduce  to 
a  formula  as 


r  mr^+t 

0=  Mt+i 


(51) 


The  dimensionless  cooling  rate  without  heat  sinks  can  be  ob¬ 
tained  by 


Q‘  =  T‘  -  1  (52) 

The  dimensionless  cooling  power  Qj  and  cooling  rate  of  the  unit 
Qo  versus  the  dimensionless  fluid  temperature  Fx  along  with  the 
COP  are  presented  in  Fig.  12.  The  cross  point  in  the  figure  is  found 
to  be  =  1.2,  which  is  a  design  point  as  the  limit  of  use  of  the 
TEC.  If  the  dimensionless  fluid  temperature  F  is  higher  than  the 
cross  point,  there  is  no  justification  for  use  of  the  TEC  although 
the  TEC  still  functions.  This  cross  point  is  defined  as  the  maximum 
dimensionless  temperature  T^  max.  There  is  also  a  minimum  point 


Fig.  12.  Dimensionless  cooling  power,  cooling  rate  of  the  unit,  and  COP  versus 
dimensionless  fluid  temperature.  This  plot  was  generated  with  Nk  =  0.3,  Nh  =  l, 
N,  =  0.5  and  ZT„2  =  1.0. 


H.  Lee/ Applied  Energy  106  (2013)  79-88 


Nh 


Fig.  13.  Optimal  (cooling  power  optimized)  dimensionless  cooling  power  and  COP 
versus  dimensionless  convection.  This  plot  was  generated  with  T/  =  0.967  and 
ZT^a  =  1.0. 


ri2h2A2  (W/K) 


Fig.  14.  Optimal  cooling  power  versus  convection  conductance  in  fluid  2  as  a 
function  of  dimensionless  convection.  This  plot  was  generated  with  T^2  =  30  °C, 
T;  =  0.967  and  ZT^2  =  1.0. 

at  =  0.83,  where  Qj  =  0,  which  is  called  the  minimum  dimen¬ 
sionless  temperature  T^  min.  Note  that  the  TEC  can  perform  effec¬ 
tive  cooling  within  a  range  from  T^ jlnta  =  0.83  to  T^  max  =  1 .2. 

Since  there  is  an  optimal  combination  of  N,  and  Nk  for  a  given 
Nh,  we  can  plot  the  optimal  dimensionless  cooling  power  Qj  opt 
and  COPopt  as  a  function  of  dimensionless  convection  Nh,  which  is 
shown  in  Fig.  13.  It  is  seen  that  both  Qjopt  and  COPopt  increase 
monotonically  with  increasing  Nh.  According  to  Eq.  (35),  the  actual 
optimal  cooling  power  Q,  opt  is  the  product  of  Qj  opl  and  r\2h2A2, 
seemingly  increasing  linearly  with  t]2h2A2.  In  practice,  there  is  a 

Table  2 

Inputs  and  results  from  the  dimensional  analysis  for  a  thermoelectric  cooler  module. 


controversial  tendency  that  Nh  may  decrease  systematically  with 
increasing  rj2h2A2.  As  a  result  of  this,  it  is  needed  to  examine  the 
variety  of  the  rj2h2A2  as  a  function  of  Nh  for  the  optimal  cooling 
power.  Fig.  14  reveals  the  intricate  relationship  between  rj2h2A2 
and  Nh  (or  r / , h , /I T )  along  with  the  optimum  actual  cooling  power 
(not  dimensionless),  which  would  lead  system  designers  to  a  vari¬ 
ety  of  possible  allocations  ( rj2h2  A2  and  Nh )  for  their  optimal  design. 
Note  that  the  analysis  so  far  is  entirely  based  on  the  dimensionless 
parameters.  Now  we  look  into  the  actual  optimal  design  with  the 
actual  values. 

For  example,  using  Figs.  13  and  14,  we  develop  an  optimal  de¬ 
sign  for  an  automobile  air  conditioner.  A  thermoelectric  cooler 
module  with  a  5-cm  x  5-cm  base  area  is  subject  to  cabin  air  at 
20  °C  in  fluid  1  and  ambient  air  at  30  °C  in  fluid  2.  We  estimate 
an  available  maximum  convection  conductance  in  fluid  1  (cabin 
air)  with  =  0.8,  hr  =  60  W/m2  K,  and  Ar  =  1000  cm2  and  also  an 
available  maximum  convection  conductance  in  fluid  2  (ambient 
air)  with  r\2  =  0.8,  h2  =  60  W/m2  K,  and  A2  =  1000  cm2,  which  gives 
ijjhjAj  =  t]2h2A2  =  4.8  W/K  and  Nh  =  1.  Note  that  >71  and  rj2  are  the 
fin  efficiencies,  and  h]  and  h2  are  the  reasonable  convection  coeffi¬ 
cients  for  the  cabin  air  cooling  and  the  ambient  air  cooling,  respec¬ 
tively.  Each  area  of  A r  and  A2  is  based  on  20  fins  (two  sides)  with  a 
fin  height  of  5  cm  on  a  5-cm  x  5-cm  base  area  of  the  module.  The 
typical  thermoelectric  material  properties  are  assumed  to  be  ocp  = 
-0in  =  220  pV/K,  pp  =  pn  =  1 .0  x  1 0  3  Cl  cm,  and  kp  =  kn=  1.4 
x  10  2  W/cm  K.  The  above  data  approximately  determines  three 
dimensionless  parameters  as  Nh  =  1,  T^  =  0.967  and  ZTrx2  =  1.0. 

As  mentioned  before,  the  present  dimensional  analysis  enables 
the  three  dimensionless  parameters  (Nh  =  l,  Tx  -  0.967  and 
ZToo2  =  1 .0)  to  determine  the  rest  optimal  parameters,  which  are 
found  to  be  Nk  =  0.3  and  N/  =  0.5  as  shown  before.  This  leads  to  a 
statement  that,  if  two  individual  fluid  temperatures  on  heat  sinks 
connected  to  a  thermoelectric  cooler  module  are  given,  an  opti¬ 
mum  design  always  exists  with  the  feasible  mechanical  constraints 
that  present  Nh.  This  optimal  design  is  indicated  approximately  at 
Point  1  in  Fig.  14.  Note  that  there  are  several  ways  to  improve  the 
optimal  power  output  by  increasing  either  rj2h2A2  or  Nh  or  both, 
which  apparently  depends  on  the  feasible  mechanical  constraints, 
whichever  is  available. 

The  inputs  and  optimum  results  at  Point  1  in  Fig.  14  are  summa¬ 
rized  in  the  first  two  columns  of  Table  2.  The  inputs  are  the  geom¬ 
etry  of  thermocouple,  the  material  properties,  two  fluid 
temperatures,  and  the  available  convection  conductance  in  fluid 
2.  The  dimensionless  results  are  converted  to  the  actual  quantities. 
The  optimal  cooling  power  is  found  to  be  54.4  W  for  the 
5  cm  x  5  cm  base  area  of  the  module.  The  cooling  power  density 
is  calculated  to  be  2.18  W/cm2.  When  iith^A^  is  limited,  simply 
increasing  tj2h2A2  invokes  decreasing  Nh,  which  is  shown  in 
Fig.  15.  This  attributes  to  the  heat  balance  that  the  hot  and  cold 
junction  temperatures  must  decrease  when  more  heat  is  extracted 


Input 

Qj  opt  (dimensionless) 

Qi  opt  (actual) 

COP i/2opt  (dimensionless) 

C0P„2opt  (actual) 

T,xl  =  20  °C,  Trxl  =  30  "C 

Nk  -  0.3 

n  =  257 

IV*  =  0.15 

n  =  128 

A  =  2  mm2,  L  =  lmm 

Nh-1 

rj^hp t,=4.8W/K 

JVh  =  l 

il,h,A,  =  4.8  W/K 

*! 2  =  0.8,  h2  =  60W/m2K, 

N,  =  0.5 

/  =  6.36  A 

Nt  =  0.25 

1  =  3.18  A 

A2  =  1000  cm2 

V!X  =  0.967 

T«,i  =  20  °C 

C  =  0.967 

Tool  =  20  ”C 

Base  area  Ab  =  5  cm  x  5  cm 

zr„,2  =  i.o 

zr„,2  =  i.o 

ZT,i2  =  1.0 

ZToc2  =  1.0 

ap  =  -a„  =  220  pV/K 

T(  =  0.930 

T,  =  8.7  “C 

TJ  =  0.949 

T,  =  14.4  "C 

pp  =  p„  =  1.0xl0  3Ocm 

T2  =  1.145 

T2  =  73.9  °C 

T2  =  1.031 

T2  =  39.4  “C 

kp  =  k„  =  1.4  x  10-2  W/cm  K 

Q5  =  0.037 

Q,  =  54.4  W 

Qj  =  0.019 

Q,  =  26.9  W 

(Z  =  3.457  x  10  3K  ’) 

COP  =0.35 

COP  =  0.35 

COP  =1.49 

COP  =1.49 

(R  =  0.01  Cl  per  thermocouple) 

Nv=  0.715 

V„  =  24.5  V 

Nv=  0.332 

Vn  =  5.7  V 

TLma*  =  1-2 

Tool,™*  “91-3  °C 

Too, max  =  3  -06 

Tooi.ma*  =  49.7°C 

T«,min  =  0-83 

Tool , min  =  -21.8  °C 

C,min  =  0-84 

Tooi,min  =  -19.2  °C 

(Cooling  power  density  Pd  =  Qj/Ab) 

Pd  =  2.18  W/cm2 

Pd  =  1.08  W/cm2 

Lee/ Applied  Energy  106  (2013)  79-88 


Nh 


Fig.  15.  Two  junction  temperatures  and  cooling  power  versus  convection  conduc¬ 
tance  in  fluid  2.  This  plot  was  generated  with,  =  4.8  W/K  T„,2  =  30  “C, 

=  0.967  and  ZJ„a  =  1.0. 

from  the  limited  input,  which  is  a  characteristic  of  the  thermoelec¬ 
tric  cooler  with  heat  sinks.  The  hot  and  cold  junction  temperatures 
are  sometimes  a  design  factor,  noting  that  the  optimal  cold  junc¬ 
tion  temperature  Tli0pt  reaches  zero  Celsius  at  Nh  =  0.4,  which 
may  cause  icing  and  reducing  the  heat  transfer. 

As  mentioned  before,  the  optimal  COP  is  sometimes  in  demand 
in  addition  to  the  optimal  cooling  power.  However,  the  real  opti¬ 
mal  COP  usually  gives  a  very  small  value  of  the  cooling  power,  al¬ 
beit  the  high  COP.  Therefore,  in  the  present  work,  the  midpoints  of 
the  optimal  N,  and  Nk  are  used  to  provide  approximately  a  half  of 
the  optimal  cooling  power  and  at  least  four  folds  of  the  cooling- 
power-optimized  COP.  This  modified  optimal  COP  is  called  a  half 
optimal  coefficient  of  performance  COP1/2opt-  These  results  are  also 
tabulated  in  the  last  two  columns  of  Table  2,  so  that  designers 
could  determine  which  optimum  is  better  depending  on  the  appli¬ 
cation.  The  optimum  cooling  power  is  usually  selected  when  the 
resources  (electrical  power  or  capacity  of  coolant)  are  abundant 
or  inexpensive  or  the  efficacy  is  not  important  as  in  microprocessor 
cooling,  while  the  half  optimum  COP  is  selected  when  the  re¬ 
sources  are  limited  or  expensive  or  the  efficacy  is  important  as  in 
automotive  air  conditioners. 

3.  Conclusions 

The  present  paper  presents  the  optimal  design  of  thermoelectric 
devices  in  conjunction  with  heat  sinks  introducing  new  dimen¬ 
sionless  parameters.  The  present  optimum  design  includes  the 
power  output  (or  cooling  power)  and  the  efficiency  (or  COP)  simul¬ 
taneously  with  respect  to  the  external  load  resistance  (or  electrical 
current)  and  the  geometry  of  thermoelements.  The  optimal  design 
provides  optimal  dimensionless  parameters  such  as  the  thermal 
conduction  ratio,  the  convection  conduction  ratio,  and  the  load 
resistance  ratio  as  well  as  the  cooling  power,  efficiency  and  high 
and  low  junction  temperatures.  The  load  resistance  ratio  (or  the 
electrical  current)  is  a  well  known  characteristic  of  optimum  de¬ 
sign.  However,  it  is  found  that  the  load  resistance  ratio  would  be 
greater  than  unity  (1.7  in  the  present  case).  This  is  a  confusing  fac¬ 
tor  in  optimum  design  with  a  TEG.  One  should  not  assume  that  the 
load  resistance  ratio  RJR  is  equal  to  unity  for  a  TEG  with  heat 
sink(s).  The  optimal  thermal  conductance  [(n)(A/I)(k)]  consists  of 
the  number  of  thermocouples,  the  geometric  ratio,  and  the  thermal 
conductivity.  It  is  important  that  there  is  an  optimum  number  of 
thermocouples  n  for  a  given  the  convection  conductance  >72^2^2 
if  the  optimal  thermal  conductance  A/L  is  constant  or  vice  versa. 
These  are  the  optimum  geometry  of  thermoelectric  devices.  The 
information  of  the  optimal  thermal  conductance  is  particularly 


important  in  design  of  microstructured  or  thin-film  thermoelectric 
devices.  Furthermore,  there  is  a  potential  to  improve  the  perfor¬ 
mance  or  to  provide  the  variety  of  the  geometry  by  reducing  the 
thermal  conductivity. 

Finally,  it  is  stated  from  the  present  dimensional  analysis  that,  if 
two  individual  fluid  temperatures  on  heat  sinks  connected  to  a  ther¬ 
moelectric  generator  or  cooler  are  given,  an  optimum  design  always 
exists  and  can  be  found  with  the  feasible  mechanical  constraints. 
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